Systematic ab initio study of the phase diagram of epitaxially strained SrTiOs 
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We use density-functional theory with the local-density approximation to study the structural 
and ferroelectric properties of SrTiOa under misfit strains. Both the antiferrodistortive (AFD) and 
ferroelectric (FE) instabilities are considered by calculating all the phases predicted by Pertsev 
et al. [Phys. Rev. B 61, R 825 (2000)] based on the phenomenological Landau theory. The 
rotation of the oxygen octahedra and the movement of the atoms are fully relaxed within the 
constraint of a fixed in-plane lattice constant. We find a rich misfit strain-induced phase transition 
sequence which is in overall agreement with the prediction by Pertsev et al. and is obtained only 
when the AFD distortion is taken into account. Nevertheless, the calculated locations of the phase 
boundaries are different from the prediction by Pertsev et al. We also find that compressive misfit 
strains induce ferroelectricity in the tetragonal low temperature phase only whilst tensile strains 
induce ferroelectricity in the orthorhombic phases only. The calculated FE polarization for both 
the tetragonal and orthorhombic phases increases monotonically with the magnitude of the strains. 
The AFD rotation angle of the oxygen octahedra in the tetragonal phase increases dramatically as 
the misfit strain goes from the tensile to compressive strain region whilst it decreases slightly in the 
orthorhombic (F04) phase. This reveals why the polarization in the epitaxially strained SrTiOs 
would be larger when the tensile strain is applied, since the AFD distortion is found to reduce the 
FE instability and even to completely suppress it in the small strain region. Finally, our analysis 
of the average polar distortion and the charge density distribution suggests that both the Ti-O and 
Sr-O layers contribute significantly to the FE polarization. 

PACS numbers: 77.22.Ej, 77.84.Dy, 77.80.-e 



I. INTRODUTION 

The strontium titanate (SrTiOa, STO) is a prototypi- 
cal perovskite compound. A mechanically free bulk STO 
will go through two soft mode-driven phase transitions as 
the temperature is lowered. First, there is the cubic-to- 
tetragonal antiferrodistortive (AFD) phase transition at 
105 K 1, 2], which results from the phonon induced insta- 
bility of the R point of the Brillouin zone (BZ) boundary 
and involves the rotation of the TiOe octahedra. This is 
then followed by the ferroelectric (FE) phase transition 
driven by the instability of the zone-center soft phonon 
mode. However, the FE transition is not observed in 
pure STO due to quantum fluctuations Therefore, it 
is regarded as an incipient ferroelectric. 

In the perovskite ferroelectrics, it is well known both 
experimentally and theoretically that the polarizations 
are strongly coupled to strain 4] , and thus the properties 
such as FE transition temperature and polarization mag- 
nitude are quite sensitive to external stress. Recently, 
using techniques such as molecular beam epitaxy and 
pulsed-laser deposition, it has become possible to grow 
high-quality perovskite films, and has attracted a great 
deal of attention. In particular, Haeni et al. recently 
reported a giant shift in FE transition temperature to 
293 K for a STO film under biaxial strains. They grew 
STO films on a (llO)-oriented DyScOs substrate, corre- 
sponding to a biaxial tensile strain of 0.8 %. The results 
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indicate that even under small strains, a nonpolar STO 
film can be dramatically transformed to a polar state. 

Several theoretical studies have provided insight into 
the structural properties and temperature dependence of 
phase transitions in an epitaxial STO film, [a, H H 01 A 
phenomenological study of STO thin films under strain 
by Pertsev et a/.|3|, based on a Landau model fitted 
to experimental data from bulk phases, predict a rich 
temperature-strain phase diagram. The effects of in- 
plane strain and an applied electric field on the dielectric 
properties of STO epitaxial thin films were calculated by 
Antons et al. 7] using density-functional theory with the 
local-density approximation (LDA). They find that the 
polarization can be tuned by the in-plane lattice param- 
eter. Recent ah initio calculations by Hashimoto et a/. 01 
confirmed that the ferroelectricity could be induced by 
applying biaxial strains and also indicated a mechanical 
origin of the ferroelectricity. 

However, in all previous first-principles calculations 
on the expitaxial STO films, the AFD instability is ne- 
glected. This is thought of as an approximation to the 
high temperature phases. Consequently, many epitaxial 
strain induced low temperature phases were unable to be 
investigated in the previous calculations. In this paper, 
we investigate the effects of the misfit strains on the struc- 
ture and ferroelectricity in the epitaxially strained STO 
with both the FE and AFD instabilities taken into ac- 
count. We use parameter-free first-principles LDA calcu- 
lations to systematically analyze all the predicted phases 
by Pertsev et al. ^ As a result, we find a rich epi- 
taxial strain induced low temperature phase transition 
sequence that was unobtainable in the previous ah initio 
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calculations. We further show the effects of the strains 
on the AFD instability itself, and hence the difference in 
the orientation of the oxygen octahedra and also the po- 
larization magnitude compared with the case when the 
AFD instability is ignored. 

The rest of this paper is organised as follows. In Sec. 
II, we describe the computational technique used and 
summarize the salient results of the previous phenomeno- 
logical analysis. In Sec. Ill, we present and discuss the 
results of our calculations for the epitaxial STO films un- 
der misfit strains. We begin with the determination of 
the theoretical stable structures of all the possible phases, 
and then proceed to relate the FE and AFD instabilities 
to the epitaxial strain and also the underlying electronic 
band structure. Finally, we give a summary in Sec. IV. 



II. THEORY AND COMPUTATIONAL 
METHOD 

Our ab initio calculations for the STO were per- 
formed using accurate full-potential projector augmented 
wave method, as implemented in the 

VASP[1J, lia package. They are based on density- 
functional theory with the LDA. PAW potentials are used 
to describe the electron-ion interaction, with the 4s, 4p, 
M and 5s states of Sr, 3d and 4s states of Ti, and 2s, 
2p of O being treated as valence states. A conjugate- 
gradient minimization scheme is used to minimize the 
Kohn-Sham energy functional, with a plane-wave cut- 
off energy of 400 eV used throughout. All calculations 
are performed with a 6x6x6 Monkhorst-Pack /c-point 
mesh. To establish minimum-energy configurations the 
total forces acting on the atoms are converged to less 
than 0.01 eV/A. Our calculations are carried out at the 
theoretical equilibrium lattice constant of 3.87 A which is 
^--^0.9 % less than the experimental value 3.905 A, under 
different strains between 2% and —2%. 



A. Structure optimization 

We begin by systematically performing optimizations 
of the five-atom unit cell of the high-temperature phases 
(FTl, HT and FOl) and the twenty-atom supercell of 
all the other possible phases involving the rotations of 
the oxygen octahedra proposed by Pcrtsev et al.^ In 
Fig. 1, we illustrate the supercell which is chosen to 
describe the AFD instability and give a schematic rep- 
resentation of the atomic displacements for the R point 
[(lll)7r/a] phonon mode. The epitaxial strain on the 
structural properties of STO is modeled by relaxing ei- 
ther the five-atom or twenty-atom cell of bulk STO under 
the constrained in-plane lattice constants, differing from 
the theoretical cubic lattice constant by fractions rang- 
ing from -2 % to 2 %. A description of these phases is 
given in Table I. We have worked out the space groups of 
all the phases which are also listed in Table I. For each 



phase, starting from a structure in which the symmetry 
is established by appropriately displacing the Ti and O 
atoms, we relax the atomic positions until the atomic 
forces fall below a small threshold of 0.01 eV/A. 

Na Sai and Vanderbilt jfj] have considered the possible 
rotations of the oxygen octahedra in the STO but only 
at zero strain. Here we consider the possible rotations of 
the oxygen octahedra in the STO under the misfit strains. 
We focus on the ground state tetragonal phase obtained 
by freezing in an AFD phonon mode at the R point of 
the Brillouin zone (BZ) boundary. This phonon mode 
corresponds to the rotation of the TiOe octahedra in the 
opposite directions from one unit cell to the next. For 
the order parameter gi, we take, e.g., the a;-axis as the 
rotation axis, and adopt a tetragonal supercell with the 
lattice vectors of length \/2a, V^a, and c along the [Oil], 
[Oil], and [200] directions, respectively. Similarly, cychc 
permutations of the lattice vectors would give rise to the 
other two rotation axes. 



B. Thermodynamic theory analysis 

Pertsev et al. used a phenomenological Landau the- 
ory to describe the mechanical substrate effect on equilib- 
rium states and phase transitions in epitaxial STO thin 
films, and developed a misfit strain-temperature phase 
diagram. Instead of dealing with thin films directly, they 
computed the structure of the bulk material strained to 
match a given cubic substrate with square surface sym- 
metry. On the basis of a Landau theory fitted to ex- 
perimental data from bulk phases, the thermodynamic 
description may be developed from the power-series ex- 
pansion of the Helmholtz free-energy density F in terms 
of temperature T, polarization P and order-parameter q, 
which corresponds to linear oxygen displacements corre- 
lated with the rotation of the oxygen octahedra around 
one of fourfold symmetry axes, and lattice strain S, de- 
fined as S = (a — ao)/a (where a is the effective lattice 
constant of the substrate and oq denotes the lattice con- 
stant of free-stranding cubic cell). By minimizing the 
free energy F{T, S*, P, q) with respect to the components 
of the polarization and structural order parameter, they 
predicted the equilibrium thermodynamic states of STO 
films under various temperature-strain conditions. The 
misfit strain-temperature phase diagram of STO is more 
complex than BaTiOa (BTO) because the coexistence of 
two coupled AFD and FE instabilities in this crystal. The 
results are listed in Table|l| Apart from the high temper- 
ature tetragonal (HT) phase, they contain purely struc- 
tural tetragonal and orthorhombic states (ST and SO), 
purely ferroelectric tetragonal and orthorhombic phases 
(FTl and FOl), and four coupled states (FT2, F02, F03 
and F04). 

In this paper, using parameter-free total-energy meth- 
ods based on density functional theory, we map out the 
equilibrium structures of the STO as a function of epi- 
taxial strain at zero temperature, and find out the effect 
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of misfit strain on the magnitude and orientation of the 
polarization. The electric polarization is calculated using 
the Berry phase method [l5L 

III. RESULTS AND DISCUSSION 

Because of the competition between ferroelectric and 
structural instabilities, the phases in the STO arc more 
complex than the BTO case0|. For clearly displaying 
the calculated energy for each phase as a function of the 
misfit strain, we separate them into two classes, which are 
separated by a large energy difference, as shown in Figs. 2 
and 5. One contains the high temperature phases, which 
do not involve the rotation of the oxygen octahedra (FTl, 
HT, and FOl). The other includes the low temperature 
phases with the coupling between the polarization and 
the rotation of the oxygen octahedra. 



A. High temperature phases 

The total energy for each of the three high energy 
phases as a function of misfit strain is shown in Fig. 2. 
For large compressive strains, the stable structure is the 
FTl phase; for large tensile strains, the FOl phase is 
prefered. At the misfit strain of around -0.1 %, there is 
a continuous transition from the polar FTl phase to the 
nonpolar HT phase. As the in-plane misfit strain is in- 
creased to 0.3 %, the STO transits from the HT phase 
to the polar FOl phase. In Fig. 3, we show the cal- 
culated polarizations for the FTl and FOl phases versus 
the misfit strain. The polarization increases dramatically 
with the magnitude of the misfit strain. For compressive 
strains larger than -0.1 %, the FTl phase has the po- 
larization along [001] whilst the FOl has zero polariza- 
tion. As the misfit strain gets close to zero, a so-called 
paraelectric gap within which all the phases exhibit zero 
polarization, appears (Fig. 3). For tensile strains larger 
than 0.3 %, the FOl phase acquires a large polarization 
along [110], whilst the polarization of the FTl phase be- 
comes zero. 

In Fig. 4, we show the calculated atomic displacements 
as a function of misfit strain. Under large compressive 
strains, the atoms in the FTl phase relax only along 
the [001] direction. As the in-plane strain decreases, a 
second-order FTl HT phase transition follows, and 
the magnitude of the atomic displacements along [001] 
gradually diminishes. While the displacements along 
[001] in the FTl phase vanish around the zero strain, 
the displacements in the [110] direction in the FOl be- 
gin to grow as we enter the tensile strain region. With 
the increasing tensile strain, the displacements in the xy 
plane in the FOl phase continue to grow monotonically. 
Near zero misfit strain, the atomic displacements change 
from along the z-direction to along an xy plane direction, 
and the change of the orientation of the polarization fol- 
lows. From these results, we can see that the in-plane 



misfit strain make the polarization to be a monotoni- 
cally increasing function of the magnitude of the strain, 
though the polarization direction switches in the para- 
electric gap. 

The above results are similar to that of the ab initio 
LDA calculations in Ref. 0- Nevertheless, there are some 
quantitative differences. For example, our estimation of 
the paraelectic gap (between -0.1 % and -1-0.3 % ) is nar- 
rower than that previously obtained (between -0.75 % 
and -1-0.54 % in Ref. 7), and this may be attributed to 
the different potentials used in the present and previous 
calculations. Furthermore, because the rotations of the 
oxygen octahedra were not included in the analysis in 
Ref. 7, the results of Ref. correspond to our results for 
the high temperature FOl, HT and FTl phases only. 



B. Low temperature phases 

Due to the interaction between the FE polarization 
and the AFD structural instability, the low temperature 
phases are more complex than the high temperature ones. 
In Fig. 5j^ we show all the low energy phases predicted 
in Ref. 6 for the misfit strains between -2 % and 2 %. 
Roughly speaking, we get nearly the same features as the 
high temperature case: the presence of two wide misfit 
strain ranges in which the STO becomes a ferroelectric. 
This is due to the coupling between the polarization and 
the misfit strain. In the films grown on the compressive 
substrates, a FE phase with polarization along the out-of- 
plane orientation ([001]) appears, while in the case of the 
tensile substrates, the polarization direction becomes in- 
plane ([110]). Nevertheless, there are some subtle differ- 
ences between the high and low temperature cases. First, 
the structures in the low temperature cases are more sta- 
ble than in the high temperature ones, because they have 
a significantly lower total energy (Figs. 2 and 5). Under 
large compressive strains, the FT2 phase has the lowest 
energy among the low temperature phases. When the 
misfit strain is decreased, the STO transits from the polar 
FT2 phase to the nonpolar ST phase and the two energy 
curves overlap in a rather wide range of tensile strains 
(Fig. 5). In contrast, for the high temperature phases, 
the energy curve of the nonpolar HT phase overlaps with 
the polar FTl phase in the tensile strains only, but with 
the polar FOl phase in both the tensile and compressive 
strains. For large tensile strains, the F04 phase is en- 
ergetically most favorable. There are some other stable 
phases occuring only in a small strain region near zero. 
The energy curves of the F02 and SO are almost the 
same, and their trends are also similar to the F04 phase 
in the small strain region, but split off from the F04 with 
increasing tensile strain (Fig. 5). Similarly, the curve of 
the F03 phase lies on the ST phase one in the compres- 
sive strain but splits off from it when grown on the tensile 
substrate. As the compressive or tensile strain increases, 
the structural symmetry is broken, thereby splitting the 
energy degeneracy. However, in the small strain region. 
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the total energies for many different phases are almost 
in overlap with each other. Nevertheless, accurate dis- 
tinction between different stable phases can be achieved 
by calculating the second derivative of the total energy 
with respect to the strain. The strain divisions obtained 
this way, between FT2 and ST, ST and F03, F03 and 
F02, F02 and F04 are -0.558 %, -0.221 %, 0.197 %, and 
0.442 %, respectively. Remarkably, this rich low tem- 
perature phase-versus-strain sequence is generally con- 
sistent with the prediction of the phenomenological the- 
ory by Pertsev et al. |^ The main difference is the precise 
location of the low temperature phase boundaries. In 
particular, we predict that the F03 phase is located in 
the strain region between -0.22 % and 0.20 % whilst it 
falls in the strain region from around -0.023 % to -0.005 
% according to Pertsev et al. &] The F02 is predicted 
to stride the zero strain in Ref. instead of the F03 
phase. Since the F02 and F03 phases have different 
polarization directions and also different order parame- 
ter q, this difference can perhaps be easily examined by 
further experiments. Of course, our calculations are for 
the zero and very low temperatures only. Consequently, 
we cannot discuss the temperature effect on the epitax- 
ial STO structures which may lift the energy degeneracy, 
and hence the misfit strain-temperature phase diagram. 

In Fig. 6, we show the calculated polarization (in- 
cluding both ionic and electronic contributions) of the 
FT2, F02, F03 and F04 phases as a function of mis- 
fit strain. Interestingly, the energy curves of all the low 
temperature phases converge to the same energy at the 
zero strain (Fig. 5). However, away from this narrow 
window near the zero strain, the lowest energy phase is 
the FT2 under the large compressive strains and the F04 
in the tensile strain region (Fig. 5). Therefore, we con- 
centrate on these two phases only below. The results are 
similar to the high-temperature phase one. The polar- 
ization is along [001] and [110] for the compressive and 
tensile strains, respectively. It exhibits a monotonic in- 
crease with the magnitude of the misfit strain. There is 
also a paraelectric gap centered at the zero strain, now 
between -0.4 % and 0.4 %. This gap is larger than the 
high temperature one. Nevertheless, the size of this gap 
should be the upper bound because we have not con- 
sidered the other FE phases in the paraelectric region. 
The polarization for the F02 and F03 phases becomes 
zero when the misfit strain drops below 0.2 % and 0.3 
%, respectively. Therefore, the overall paraelectric gap 
for the low temperature phases is between -0.4 % and 
0.2 %. Interestingly, a comparison of the polarization- 
versus-strain curves for the FTl phase in Fig. 3 and the 
FT2 phase in Fig. 6 reveals that the AFD instability 
reduces the polarization in the large compressive strain 
region, and completely suppresses it in the small strain 
region, thus enlarging the paraelectric gap. Note that the 
atomic structures of the FTl and FT2 phases are almost 
the same except that the FT2 phase can accommodate 
the AFD rotations of the oxygen octahedra whilst the 
FTl phase cannot. Indeed, the AFD distortion as char- 



acterized by the rotation angle of the oxygen octahedra 
in the FT2 is strong in the compressive strain region, as 
shown in Fig. 8 below. 

We now further quantify the strength of the FE and 
AFD distortions. In the FE case, let us define the average 
polar distortion along the z-direction, p^, in the 1x1 unit 
cell of each atomic layer. For the Ti02 layers. 



Pi = Ai{Ti) 
and for the SrO layers, 



A,(0i)+A,(02) 



where are the atomic displacements in the i-direction 
relative to the undistorted symmetric 1x1 perovskite 
structure. The computed results describing the strength 
of the FE distortion in each atomic layer are ploted in 
Fig. 7. Note that for the ST and SO phases, the average 
polar distortion is zero for all the strains, and therefore, 
is not shown in Fig. 7. 

It is clear from Fig. 7 that in the compressive strain 
region, the average polar distortion in the FT2 phase is 
larger in the Ti02 layers than in the SrO layers. Inter- 
estingly, in the F04 phase, it is the opposite, i.e., the 
average distortion in the tensile strain region is larger in 
the SrO layer than in the Ti02 layer. Around the zero 
strain, the average distortion is significant only in the 
Ti02 layer in the F03 and F02 phases. This indicates 
that in the vicinity of the zero strain, only the F02 and 
F03 phases have a significant polarization, thereby giv- 
ing a microscopic explanation to the results of the previ- 
ous thermodynamical analysis ^3] (Table I). Note also that 
the average polar distortion along the in-plane direction 
in the tensile strain region in the F03 and F04 phases 
is larger than along the z-direction in the compressive 
strain region in the FT2 phase. This suggests that the 
ionic contribution to the polarization in the F04 phase 
in the tensile strain region would be larger than in the 
FT2 phase in the compressive strain region, and explains 
the calculated polarizations displayed in Fig. 6. 

The strength of the AFD distortion is characterized 
by the Ti04 rotation angle 9, in each Ti02 layer. In 
Fig. 8, we present the calculated rotation angle of the 
oxygen octahedra as a function of misfit strain for the 
low-temperature phases due to the AFD instability. The 
calculated tilting angle of the oxygen octahedra at zero 
strain (6°) is overestimated in our calculations, as also in 
Ref. Q, compared to the experimental value of 2°[l8ll. 
This discrepancy between the ab initio calculations and 
the experiments may be caused by the LDA error or by 
the quantum fluctuations 1^ which should reduce the 
rotations of the oxygen octahedra, or by both. The re- 
sults nonetheless provide us some qualitative insight into 
how oxygen octahedra rotate under the applied strains. 
Fig. 5 shows that under compressive strains, the FT2, 
ST and F03 are the low energy phases and the rotation 
angle of the TiOe octahedra along [001] increases with 
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the magnitude of the misfit strain. In contrast, under 
the tensile strains, the F04, F02 and SO become the 
low energy phases (Fig. 5) and the rotation axis changes 
to lie in the xy plane. In the F02 and SO phases, the ro- 
tation angle of the TiOg octahedra again increases with 
the magnitude of misfit strain (Fig. 8). Interestingly, 
the rotation of the oxygen octahedra in the F04 phase 
is nearly independent of appHed strain (Fig. 8). Note 
that there are two equivalent rotation axes along [100] 
and [010], respectively, for the F04 phase. Therefore the 
linear displacements of the oxygen atoms due to the rota- 
tions and the rotation angles are equal for the two axes. 
In Fig. 8, we thus show only the rotation angle for either 
[100] or [010]. Interestingly, Figs. 6 and 7 show that the 
FE instability as manifested in the electric polarization 
and the average polar distortion, in both the FT2 phase 
in the large compressive strain region and the F04 phase 
in the large tensile strain one, grow as the magnitude of 
the strains is increased, though it is stronger in the F04 
phase than in the FT2 phase. In contrast, the behavior 
of the AFD instability as manifested in the rotation of 
the oxygen octahedra, is different. Fig. 8 shows that the 
tilting angle of the oxygen octohedra in the FT2 (also 
ST and F03) phase increases dramatically as the misfit 
strain goes from the strong tensile region through zero 
strain to the strong compressive region, while it remains 
nearly as a relatively small constant for the F04 phase. 
This explains why the calculated polarization is larger 
for the F04 phase under the tensile strain than the FT2 
phase when the AFD instability which tends to suppress 
the FE instability, is taken into account. 



C. Charge density and polarization 

In order to better understand the chemical bonding 
and also the origin of the electric polarization in the 
strain induced FE phases, we examine the calculated 
charge density distributions. Figure 9 shows the contour 
plots of the charge density differences between calculated 
charge densities and those obtained by a superposition of 
the free atomic charge densities for the FT2 phase under 
compressive strain of -1.6 % and the F04 phase at tensile 
strain of 1.6 %. It can be seen from the calculated charge 
density difference distribution on the (001) SrO planes in 
Fig. 9 that in the both structures, the charge is depleted 
in the inner region of the Sr sites, whilst, in contrast, 
there is considerable charge buildup on the O sites and 
also in the outer region of the Sr sites. This indicates that 
the Sr-0 bonding is predominantly ionic. The ionic Sr-0 
bonding can be further seen from the fact that the charge 
distrubtion at the Sr sites are nearly spherical. This is 
because the Sr 5s valence electrons have moved to the 
neighboring O atoms which have high electron afRnitiy, 
with the spherical symmetry ionic cores left behind. On 
the other hand, the bonding between the Ti and O atoms 
are strongly covalent, as can be seen from the calculated 
charge density difference distributions on the (001) and 



(100) TiO planes in Fig. 9. There is pronounced charge 
buildup in the vinicity of the Ti-0 bond center whilst 
there is significant charge depletion on both the Ti and 
O sites. 

Figure 9 also shows that the charge density distribu- 
tion on both the (001) SrO and TiO planes in the FT2 
phase has a fourfold symmetry at the O and Ti sites, re- 
spectively. As a result, there is no electric polarization 
in the x-y plane. In contrast, on the (001) TiO plane 
in the F04 phase [Fig. 9(d)] the charge buildup is much 
stronger in the central region of the two right Ti-0 bonds 
than in the left two ones, and it leans towards the verti- 
cal down direction in all the four Ti-0 bond regions, thus 
resulting in an electric polarization along the [110] direc- 
tion (vertical up direction). Clearly, the charge density 
distribution on the (001) SrO plane in the F04 phase also 
lean downwards, suggesting that the SrO layers also con- 
tribute to the polarization, being in conflict with common 
wisdom. On the (100) TiO planes [Fig. 9(e-f)],the strong 
charge buildup in the FT2 phase clearly lean downwards 
and is symmetric with respect to the vertical line through 
the center O atom in the FT2 phase, whilst, in contrast, 
it leans towards left and is symmetric with respect to 
the horizontal line through the center O site in the F04 
phase. As a result, the FT2 phase has a polarization 
along [001] whilst the F04 phase exhibit a polarization 
along [110]. 



IV. SUMMARY 

In this paper, we have performed ab initio LDA cal- 
culations for the STO under various misfit strains. A 
unique part of this work is that we have considered both 
the high-temperature and low-ternperature phases of the 
STO proposed by Pertsev et al.^ such that both the 
AFD and FE structural instabilities have been investi- 
gated, unlike previous ab initio studies in which the AFD 
rotations of the oxygen octahedra have been neglected. 
As a result, we find a rich misfit strain-induced low tem- 
perature phase transition sequence which is in overall 
agreement with the prediction by Pertsev et al. Nev- 
ertheless, the calculated locations of the phase bound- 
aries are rather different from the prediction by Pertsev 
et al. In particular, we predict that the F03 phase is 
located in the strain region between -0.22 % and 0.20 
% whilst the F02 would stride the zero strain accord- 
ing to Pertsev et al. Since the F02 and F03 phases 
have different polarization directions and also different 
order parameter q, this difference can perhaps be easily 
examined by the future experiments. By analyzing the 
FE atomic displacements and the AFD rotations of the 
oxygen octahedra, we provide a comprehensive micro- 
scopic insight into the effects of both the FE and AFD 
instabilities on the ferroelectric properties of STO under 
the epitaxial strains. We also find that at low tempera- 
tures, compressive strains would induce FE polarization 
in the FT2 phase only whilst tensile strains induce po- 
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larization in the orthorhombic phases only. The calcu- 
lated ferroelectric polarization for both the FE tetrag- 
onal and orthorhombic phases increases monotonically 
with the magnitude of the strains, with the polarization 
in the tensile strain region being larger than in the com- 
pressive region. The AFD rotation angle of the oxygen 
octahedra in the FT2 phase increases dramatically as the 
misfit strain goes from the tensile to compressive strain 
region whilst it remains nearly unchanged in the F04 
phase. This reveals why the polarization in the epitaxi- 
ally strained STO would be larger when the tensile strain 
is applied, since we find that the AFD distortion would 
reduce the FE instability and even completely suppress 
it in the small strain region. Finally, by analysing the so- 



called average polar FE distortion and the charge density 
distribution, we find that both the Ti-0 and Sr-0 layers 
contribute significantly to the FE polarization. 
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Figure captions 

Fig. 1 (Color onhne) A [001] projected section of the 
SrTiOa lattice illustrating the linear displacements of 
the oxygen ions associated with an antiferrodistortive 
phonon mode at the R point [(lll)7r/a] of the Bril- 
louin zone boundary (upper panel). The thick square 
frame indicates the supercell used in our calculations. 
The lower panel is a schematic representation of the ro- 
tations of the oxygen octahedra in the soft phonon-mode. 

Fig. 2 (Color online) Calculated total energy of the 
high temperature epitaxial SrTiOa phases versus misfit 
strain (see Table for the symbols of the phases). The 
vertical dotted lines denote the boundaries between the 
paraelectric and ferroelectric phases (see Fig. 3). 

Fig. 3 Calculated polarization of the FTl and FOl 
phases versus misfit strain. The polarizations for the 
FTl and FOl phases are along [DDI] and [110], respec- 
tively. The vertical dotted lines denote the boundaries 
between the paraelectric and ferroelectric phases. 

Fig. 4 (Color online) Calculated atomic displacements 
A in the FTl phase (compressive strains only) and the 
FOl phase (tensile strains only). The subscripts label 
the directions of the displacements. 

Fig. 5 (Color online) Calculated total energy of the 
low temperature epitaxial SrTiOa phases versus misfit 
strain. The vertical dotted lines denote the boundaries 
between the paraelectric and ferroelectric phases (see 
Fig. 6). 

Fig. 6 (Color online) Calculated polarization of the 
FT2, F02, F03 and F04 phases versus different strain. 
The polarization for the FT2 phase is along [001]; the 
polarization for the F03 and F04 phases is along [110]; 
the polarization for the F02 phase is along [100]. The 
vertical dotted lines denote the boundaries between the 
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TABLE L The epitaxial SrTiOa phases predicted by Pertsev et al. Nonzero components of the polarization P and the order 
parameter q are the characteristics of the different phases in the epitaxial SrTiOs films grown on the cubic substrates. Phases 
are labeled by the lattice symmetry (T = tetragonal, O = orthorhombic) and also the type of the instability (S = structural, 
F = ferroelectric, H = high temperature). Also listed are the space groups (Gs) for all the phases found in the present work. 



Phase 


HT 


ST 


SO 


FTl 


FT2 


FOl 


F02 


F03 


F04 


P 








Pi 




\PA = \P2\ 


Pi (or P2) 


|Pl| = |P2| 


|Pi! = |P2! 


q 




93 


qi (or (72) 








12 (or gi) 


93 


kii = k2| 


Gs 


P4/mmm 


I4/mcm 


Fmmm 


P4mm 


14cm 


Amm2 


Fmm2 


lma2 


lma2 



paraelectric and ferroelectric phases. 

Fig. 7 (Color online) The average layer ferroeletric 
polar distortion (see text), pi, versus misfit strain, in 
percentage of the lattice parameter. Solid and dashed 
curves denote the distortions in the Ti02 and SrO 
layers, respectively. The subscripts label the directions 
of the distortions. Note that for the ST and SO phases, 
the average polar distortion is zero for all the strains. 

Fig. 8 (Color online) Calculated rotation angle of the 
oxygen octahedra for all the low temperature epitaxial 
SrTiOa phases versus misfit strain. 

Fig. 9 (Color online) Contour plots of the charge den- 
sity differences between the calculated charge densities 
and those obtained by a superposition of the free atomic 
charge densities in different projected planes for the FT2 
phase under compressive strain of -1.6 % and the F04 
phase at tensile strain of 1.6 %. (a) and (b) are for the 
(001) SrOa planes, (c) and (d) are for the (001) TiO 
planes, (e) and (f) are for the (100) TiO planes. The dark 
(blue) curves denote the positive contour levels with con- 
tour step of 0.025 e/A'^and the hght (red) dashed curves 
denote the negative ones with contour step of 0.084 e/K^. 
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